"0","# calculate richness and diversity"
"0","S.obs <- specnumber(pollen.6000)"
"0","H <- diversity(pollen.6000, ""shannon"")"
"0","Simp <- diversity(pollen.6000, 'simpson')"
"0","H_even <- H/log(S.obs)"
